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Abstract 

We report evidence for the decay modes X(3872) — > jJ/tp and X(3872) — ► 7r + 7r _ 7r° J/ip. In 
the latter, the 7r + 7r~7r° invariant mass distribution has a strong peak between 750 MeV and the 
kinematic limit of 775 MeV, suggesting that the process is dominated by the sub-threshold decay 
X — > uJ/ip. These results establish the charge-conjugation parity of the X(3872) as C = +1. The 
results are based on a study of X(3872) mesons produced via exclusive — > K T X (3872) decays 
in a 256 fb _1 data sample collected at the T(4S I ) resonance in the Belle detector at the KEKB 
collider. The X(3872) — > jJ/ip signal and n + ir~ it J / 'ip signal each has a statistical significance 
that is greater than 4cr. 

PACS numbers: 14.40. Gx, 12.39.Mk, 13.20.He 
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The X(3872) was discovered by Belle as a narrow 7i + 7i^J/if> mass peak in exclusive B — > 
K~ ix + it~ J / ip decays 0,0]. The observed mass and the narrow width are not compatible with 
expectations for any of the as-yet unobserved charmonium states Moreover, the tt + tt~ 
invariant mass distribution peaks near the upper kinematic limit of M(tt + 7t~) = 775 MeV, 
as expected for dipions that originate from p — > 7r + 7r - decays. Charmonium decays to 
pJ/if> final states violate isospin and are expected to be suppressed. The X(3872) and its 
above-listed properties were confirmed by other experiments 0]. 

The X[3872) mass (3871.9±0.5 MeV @) is within errors of the D°D°* threshold (3871. 3± 
1.0 MeV p); the difference is 0.6 ± 1.1 MeV. This has led to speculation that the X might 
be a molecule-like D°D°* bound state E|. According to Ref. Q, the preferred quantum 
numbers for such a bound state would be either J PC = 0~ + or 1 ++ . The decay of a 
C = +1 state to 7t + 7t~J/iIj would proceed via an isovector p°J/if> intermediate state and 
produce a it + tt~ mass spectrum that is concentrated at high masses, as is observed. In 
this meson-meson bound state interpretation, the close proximity of the X mass to D°D°* 
threshold compared to the D + D*~-D°D°* mass splitting of 8.1 MeV produces a strong 
isospin violation. 

Swanson proposed a dynamical model for the X(3872) as a D°D°* hadronic resonance 
In this model, J PC = 1 ++ is strongly favored and the wave function has, in addition to 
D°D *, an appreciable admixture of uuJ/if) plus a small pJ/ip component. The latter pro- 
duces the 7i + n~J/ip decays that have been observed; the former gives rise to n + n^n J/ip 
decays via a virtual u that is enhanced because of the large uJ/ip component to the wave- 
function. Swanson's model predicts that X(3872) — > it + J / ip decays should occur at 
about half the rate for -k + -k^ J/ip and with a tt + 7t~7t° invariant mass spectrum that peaks 
near the upper kinematic boundary of 775 MeV (7.5 MeV below the u peak). In a subse- 
quent paper jlo| . Swanson emphasized the importance of searching for the radiative decay 
X(3872) — > yJ/if). The observation of this decay mode would unambiguously establish the 
charge-conjugation parity of the X(3872) as C — +1. 

In this Letter we report on a search for B — > KjJ/if> and K^-k^-k^J /ip decays in a 
275 million BB event sample collected in the Belle detector at the KEKB energy-asymmetric 
e + e~ collider. The data were accumulated at a center-of-mass system (cms) energy of 
\fs = 10.58 GeV, corresponding to the mass of the T(4S) resonance. KEKB is described in 
detail in ref. fllj ] . 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector, a 50-layer cylindrical drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect K L mesons and to identify 
muons (KLM). The detector is described in detail elsewhere [l2T |. 

For the X(3872) — > •yJ/ip study we select events that contain a J/V, either a charged 
or neutral kaon and a 7 with laboratory-frame energy greater than 40 MeV. We use the 
J /if) and kaon selection criteria described in refs. £0] and The 7 candidate is rejected 
if, when combined with any other photon in the event, fits a 7r° — > 77 hypothesis with 
X 2 < 4. To reduce the level of e + e~ — > qq (q = u, d, s or c-quark) continuum events in the 
sample, we also require R2 < 0.4, where R2 is the normalized Fox- Wolfram moment (l4| . 
and I cos#g| < 0.8, where 8b is the polar angle of the 5-meson direction in the cms. 

Candidate B — > K^J '/if) mesons are identified by the energy difference AE = E B ms —E^ m 
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where E^ m is the 



and the beam-energy constrained mass M^ c = \J(E^ m ) 2 — (p 
cms beam energy, and E c ™ and p c ™ are the cms energy and momentum of the K^J/ip 
combination. We select events with M^ c > 5.20 GeV and \AE\ < 0.2 GeV and among these 
define a signal region 5.2745 GeV < M bc < 5.2855 GeV and \AE\ < 0.035 GeV, which 
correspond to ~ ±2cr from the central values for both variables. 

There is prominent peak near 3510 MeV in the M^J/ip) distribution [l5|, due to B — » 
Kxd] Xci — ► jJ/tp- This is used as a calibration reaction to determine Mbc AE and 
Mi^fJ/tp) resolutions and shifts in the nominal peak positions. 
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FIG. 1: The M^ c distributions for 64-MeV-wide M(jJ/tJ;) mass bins in the X(3872) mass region. 
The curves are the results of fits described in the text. 

Figures ^a) through (f) show the M^ c distributions for six 64- MeV- wide M^J/ip) mass 
bins in the V(3872) region. The 64 MeV bin- width corresponds to ±2.5cr, where a is the 
M^J/ip) mass resolution determined by scaling the measured width of the x& peak to the 
V(3872) region using a MC-determined factor. The upper right-hand panel of Fig. ^ cen- 
tered on 3872 MeV, is the only bin with a significant B meson signal. The AE distributions 
(not shown) also have a significant B meson signal in the 3872 MeV bin, and only there. 

The curves in the figures are the results of binned likelihood fits that are applied simul- 
taneously to the Mb c and AE distributions. The fit uses single Gaussians for the Mb c and 
AE signals with widths fixed at values determined from the x& iJ/ip event sample. The 
areas of the Mb c and AE signal Gaussians are constrained to be equal. For the backgrounds, 
an ARGUS function is used for Mb c and a first-order polynomial for AE. 

The results of the fits are plotted in Fig. El The signal yield in the 3872 MeV bin 
is 13.6 ± 4.4 events. We estimate the statistical significance of the excess signal in the 
M = 3872 MeV bin from J— 2 ln[£(rabkg)/£(n max )], where C(n) is the likelihood value for a 
signal of n events, nbk g is the background level in the X(3872) bin determined from a linear 
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fit to the event yields in the other bins (2.6 ± 0.6 events) inflated by la, and n max is the 
best-fit signal yield. The significance is 4.0a. 
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FIG. 2: The signal yields from fits to the M bc and AE distributions in bins of M(3872). The 
curves are results of fits described in the text. 



Using a MC-determined efficiency of (19 ± 1)%, we determine the product branching 
fraction 

B(B -> KX)B{X -> 7 J/V>) = (1.8 ± 0.6 ± 0.1) x 10" 6 , (1) 

where the first error is statistical and the second systematic. The latter includes effects of 
differences between data and MC resolutions for M bc and AE (3%), tracking efficiency (3%), 
charged kaon identification (1%), K s detection (3%) and 7 detection efficiency (2%) added 
in quadrature. We assume that charged and neutral B mesons have the same branching 
fraction to KX(3&72). The product branching fraction of Eqn. ^ corresponds to the partial 
width ratio 

T(X -> 7 J/V>)/r(X -> 7T+7T- J/ij) = 0.14 ± 0.05. (2) 

For the X — > 7r + 7r~7r° J/ip search we select events with a J/ip, a kaon, and a 7r + 7r~7r° 
combination. The J/ip, kaon, R2 and 9b requirements are the same as for the jJ/ip search 
described above. For charged pions, we use any track that is not positively identified as 
a lepton. Neutral pions are reconstructed from pairs of gamma rays detected in the Csl 
calorimeter that fit a 7r° — > 77 hypothesis with a y 2 < 6. We require the ir° cms momentum 
to be greater than 180 MeV. We eliminate events of the type B — > Ktt°iIj(2S); '0(25') — > 
7i + 7i^J/ip by vetoing events where M(7t + ti^ J/ip) is within 3a of m^s)- 

In the case of multiple entries corresponding to different tt° — > 77 assignments, we select 
the 77 combination with the best x 2 value for the 7r° — > 77 hypothesis. For multiple entries 
caused by different charged particle assignments we select the track with the smallest impact 
parameter relative to the J/ip — > £ + £~ vertex. The signal region is defined as 5.2725 GeV < 
M bc < 5.2875 GeV and \AE\ < 0.035 GeV. To select possible V(3872) -> tt+tt-tt J/^ events 
we require |AM| < 0.0165 GeV, where AM = M(3tt J/ip) - 3.872 GeV; this corresponds to 
a ±3(7 requirement. 

Figures EJa) through (1) show the M bc distributions for events in the AE and X — > 
7T + 7T~7T°J/ip signal regions for 25 MeV-wide tt + tt~7t° invariant mass bins. There is a distinct 
B meson signal for the M(7r + 7r~7r°) > 750 MeV bin (Fig. E(l)), and no evident signals for 
any of the other 3ir mass bins. The AE distributions (not shown) exhibit a similar pattern. 

The curves in the figures are the results of binned likelihood fits that are applied simul- 
taneously to the M bc and AE distributions. The fit uses single Gaussians for the M bc and 
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FIG. 3: Mb c distributions for B —* K 7r + 7r n°J/ip candidates in the AE and X — > tt + it ir°J/ip 
signal regions for 25 MeV-wide 7r + 7r~7r invariant mass bins. 



AE signals with widths fixed at their MC-determined values. The areas of the M^ c and AE 
signal Gaussians are constrained to be equal. For the backgrounds, an ARGUS function 
is used for M bc and a second-order polynomial for AE. 
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FIG. 4: The .B-meson signal yield in bins of 3n invariant mass determined from fits to the M^-AE 
distributions. 

Figure |U shows the fitted 5-meson signal yields vs M(7r + 7r~7r°). All of the fitted yields 
are consistent with zero except for the M(7r + 7r~7r°) > 750 MeV bin, where the fit gives 
12.4 ±4.1 events. 

As possible backgrounds to the observed signal, we considered feed-across from B — > 
KujJ /ip decays and non-resonant B~ — ► K~7T + 7T~7T°J/ijj decays. 
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Belle has reported the observation of a near-threshold uJ/ip mass enhancement in B — > 
Km J /if) decays |l3j . Although the peak of the uj — > 7r + 7r~7r° resonance at 782.5 MeV is 
7.5 MeV above the maximum possible 3ir invariant mass value for X — > 7t + 7t~7t°J/i/) decays, 
there is some overlap in the tails of the kinematically allowed regions for the two processes 
that might result in some events from one signal feeding into the other. From the fit results 
to the u)J/if> mass spectrum for B — > KujJ/if) of ref. k| we determine that this source 



contributes 0.75±0.14 events to the overlap region. As an independent check, we determined 
the A(3872) — > 7r + 7r~7r° J/?/> signal yield with a more severe restriction on M(7r + 7r _ 7r° J /if)), 
namely Mx — 3cr < M (-k + n~ -k° J / if)) < Mx ± lcr, that has no overlap with the uj band. The 
X — ► 7T + 7r _ 7r°J/^ signal yield in the truncated region is 10.6 ± 3.6 events. For a Gaussian 
signal distribution with no feed-across background, we expect the truncation of the signal 
region to reduce the signal by 2.1 events (16%); the observed reduction of 1.8 events is 
consistent with a feed-across level that is less than one event. 

To study possible contributions to the A(3872) — ► n + n~ -k° J / if) signal from non-resonant 
B~ — ► K~ 7t + tt~ 7i° J/ ip decays, we defined two side-regions in the [AM, M(3tc)} plane. The 
"upper" side-band is centered at AM = 0.033 GeV (just above the A(3872) signal region) 
and is ±0.0165 GeV wide, with 660 MeV < M(3tt) < 760 MeV. The "lower" sideband is 
±0.0165 GeV wide centered at AM = -0.033 GeV with M(3vr) > 640 MeV. Each sideband 
has an area that is four times that of the X — > 7r + 7r~7r° J /if) signal bin. 

There is no evidence for significant signal yields in either sideband. Fits to the M^ C -AE 
distributions give 4.3 ±6.2 and 6.4 ±5.6 signal events for the upper and lower bands, respec- 
tively. Normalizing by area, this gives an estimate of a possible non-resonant background 
level in the X — > 7i + 7i^7i°J/if) signal bin of 1.3 ± 1.0 events. 

The combined estimate for feed-across and non-resonant backgrounds is 2.1 ± 1.0. We in- 
flate this by ±lcr and compute the significance with ribkg = 3.1 events and n max = 12.4 events. 
The resulting statistical significance is 4.3a. 

To determine the branching fraction for B — > -k + n~ -k° J / if) we attribute all of the signal 
events with M(7r + 7r~7r°) > 750 MeV to X — > n + n~ -k° J / if) decay. We compute the ratio of 
7T + 7T~7T J/tp and tt + 7t~J/?I) branching fractions by comparing this to the number of X — > 
7T + 7T~J/?f) events in the same data sample, corrected by MC-determined relative detection 
efficiencies. The ratio of branching fractions is 

B(X ->• 7r + 7r-7r°JM) . , . , 

^_ T _^ =1 .0 ± 0.4(stat) ± 0.3(syst), (3) 

where the systematic error reflects the uncertainty in the relative acceptance, the level 
of possible feed-across and nonresonant background, and possible event loss due to the 
M(37r) > 750 MeV requirement, all added in quadrature. 

In summary, we report strong evidence for the decays X(3872) — > jJ/if> and X(3872) — > 
7i + 7i^7i J/ip. These are the first measurements of an X(3872) decay mode other than 
7r + 7r~ J /if). The X(3872) — > jJ/if) result unambiguously establishes the charge conjuga- 
tion parity of the X as C = ±1. This, in turn, would mean that the tt + tt~ system in 
X — > 7T + 7r _ J /if) decay is from the decay of a p meson, as is suggested by the observed 7r + 7r~ 
invariant mass distribution 

In the X — > 3irJ/if> signal, the 7r + 7r~7r° invariant masses are clustered above 750 MeV, 
near the upper kinematic boundary. This is suggestive of a sub-threshold decay to a virtual 
uj plus a J /if), and is further evidence for C = ±1. This decay mode was predicted by 
Swanson to occur at about the observed strength based^n a model where the V(3872) is 
considered to be primarily a D°D°* hadronic resonance 
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The large isospin violation implied by the near equality of the pJ/ip and LoJ/tj) decay 
widths is difficult to accomodate in a cc charmonium interpretation of the AY3872), but a 
natural consequence of the meson-meson bound state model point of view 0, @| • 
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